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Abstract: We propose a new equation of state for the Dark Energy component of the Universe. 
It is modeled on the equation of state p = w{p — p^) which can describe a liquid, for example 
water. We show that its energy density naturally decomposes into a component that behaves as 
a cosmological constant and one whose energy density scales as a~'^'^'^^'^\ and fit the parameters 
specifying the equation of state to the new SNIa data, as well as WMAP and 2dF data. We consider 
both the case where the dark fluid is smooth {i.e only the CDM component clusters gravitationally) 
as well as the case where the dark fluid also clusters. We find that for both cases, reasonable values 
of the parameters can be found that give our model the same as that of ACDM. Furthermore, 
in the case where our dark fluid clusters, allowing a blue tilt to the power spectrum allows us to fit 
all the requisite data with the dark fluid being the dominant component of the dark matter of the 
Universe. A remarkable feature of the model is that we can do all this with a microphysical w > {). 
We also display a field theoretic model that yields this equation of state. 
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The nature of the dark energy component of the Universe [1, 2, 3] remains one of the deepest 
mysteries of cosmology. There is certainly no lack of canditates: a remnant cosmological constant, 
quintessence[4], k-essence[5], phantom energy[6]. Modifications of the Friedmann equation such 
as Cardassian expansion[9] as well as what might be derived from brane cosmology[10] have also 
been used to try to explain the acceleration of the Universe. 

In this work, we offer a new candidate for the dark energy: Wet Dark Fluid (WDF). This model 
is in the spirit of the generalized Chaplygin gas (GCG) [11], where a physically motivated equation 
of state is offered with properties relevant for the dark energy problem. Here the motivation stems 
from an empirical equation of state proposed by Tait [7, 8] in 1888 to treat water and aqueous 
solutions. The equation of state for WDF is very simple, 

PWDF = w{pwBF - P*), (1) 

and is motivated by the fact that this is a good approximation for many fluids, including water, in 
which the internal attraction of the molecules makes negative pressures possible. One of the virtues 
of this model is that the square of the sound speed, c^, which depends on dp/ dp, can be positive 
(as opposed to the case of phantom energy, say), while still giving rise to cosmic acceleration in 
the current epoch. 

In real fluids negative pressures eventually lead to a breakdown of Eq. (1) due to cavitation 
[12], but for now we simply treat Eq. (1) as a phemenological equation[13]. We will show that this 
model can be made consistent with the most recent SNIa data[14], the WMAP results[15] as well 
as the constraints coming from measurements of the matter power spectrum[16]. The parameters 
w and are taken to be positive and we restrict ourselves to < ui < 1. Note that if Cg denotes 
the adiabatic sound speed in WDF, then w = c^. 

To find how the WDF energy density scales with the scale factor a, we use the energy conser- 
vation equation together with the equation of state in Eq. (1): 

PWDF + 3iJ(pwDF + PWDf) = 

^Pwu. = j^P. + DC-^f^^--\ 
1 + w a 

where -D is a constant of integration and ao is the scale factor today; we will set ao = 1 from now 
on. 

WDF naturally includes two components: a piece that behaves as a cosmological constant as 
well as a piece that redshifts as a standard fluid with an equation of state p = wp. We can show 
that ifv/e take D > 0, this fluid will never violate the strong energy condition p + p>0: 



PWDF + PWDF = (1 + w)PWDF - Wp^ 

= £,(l + ^)(^)3(l+-)>0 



a 

It is tempting to try to use the second component as the dark matter, thus unifying the two dark 
components. There is the potential problem that since the sound speed of the second component is 
non-zero, this would give rise to a pressure gradient term in the equation for linear fluctuations. At 
least in the case of the Chaplygin gas, this has been shown to slow down the growth of fluctuations 
to a level that would be inconsistent with measurements of the power spectrum[17] as well as with 
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the CMB fluctuations (although see [18, 19] for possible ways out of this predicament). We will 
start off by assuming that WDF does not cluster gravitationally and that the formation of structure 
is driven, as in the standard ACDM model, by the clustering of a cold dark matter component 
(this was done in ref.[20] for the GCG model). Later on in the paper, we will examine whether 
and to what extent fluctuations in the matter power spectrum (MPS) are suppressed if we aUow 
fluctuations in the WDF fluid and whether adjusting the tilt of the power spectrum could reduce 
such a suppresion. 

The Friedmann and acceleration equations for the WDF/CDM system in a spatially flat FRW 
universe are: 



a 2 



- = — ^ [ficDM — 2 (riwDE + S^wdm) (4) 



+ 3(1 + w)nwDM a-^(^+"')] , 

where WDE, WDM stand for wet dark energy/matter respectively and 

SttGn w SttGn „ 

^™ = ^TT^^*' ^^^^ = -jHf^' 

with riwDE + ^^WDM + ^^CDM = 1- We can also rewrite the equation of state in the more traditional 
formpwDF = Wea{z)pwDF where 

. . _ W OWDM (1 + Z)^(^+^^ - OwDE ... 
""•^^^^^ J^WDM (1 + (^+-) + J^WDE ■ ^ ^ 

We see that Wes interpolates between w at early times and —1 at late times. Note that this is of a 
similar form to that used in [21], except that some of the parameters in their Wes that they take to 
be positive are in fact negative in our version. 

There are some obvious constraints our model must satisfy, namely that the WDM component 
of WDF should not dominate the energy density of the Universe at nucleosynthesis. Comparing 
the energy density of WDM to that of radiation, we find 

PyyPM _ ^WDM ^ ~)3m>-1 
Prad ^rad 

Since the CDM and WDE components can be neglected at nucleosynthesis, we find that we can 
obtain a relation between w and Owdm by making use of [22] 

Prad = ^9*(.T)T\ (8) 

where (/*(T ~ 1 MeV) = gf'^^'^'^'^'^ + Ag^,. The contribution of standard model particles in thermal 
equilibrium at T ~ 1 MeV is ^standard ^ 10.75. Standard BBN hght element abundances lead to 
the bound ^^(r ~ 1 MeV) < 12.50, which in turn leads to: 

Ag, = + ^j^^j3«,-l^standa.d < j 75 ^ ^Wd ^ ^2.36 X 10-^0)^—^ (9) 

^ ^rad ^ 'rad 

Using the value h = 0.7 for the Hubble parameter we find that getting the correct light element 
abundances imposes the following constraint: In f^wDM < H O — 66.5 w. 
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As mentioned above, we considered three 
distinct experiments; 1) CMB measurements 
from WMAP[15], 2) matter power spectrum 
from the 2dF survey[16] and 3) type la su- 
pemovae observations[14]. We used CMB- 
Fast[24] to compute the scalar power spec- 
trum and the transfer functions for our equa- 
tion of state. The transfer functions were then 
used to obtain the matter power spectrum. 

To fit to the supemovae data we integrated 
di^z) = (l + z) Jq dz' / H{z') numerically to 
obtain the red-shift/distance modulus relation 
for our model. We then simultaneously fit all 
three experiments by constructing the com- 
bined x^- Although the full model has a num- 
ber of parameters, including many nuisance 

parameters, in practice the CMB data places strong enough constraints on the domain of J^wdm 
so that we can fix most of the parameters to their best fit value from ACDM, leaving us free to 
consider w and JIwdm- We calculated the combined on a grid of w-ilwDM values to identify 
the confidence regions and found that the fit is primarily constrained by the CMB and the matter 




Figure 1: The contours in the w-fiwDM plane. 
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ACDM 






All. 291 


Concordance model 


Set 1 


0.316228 


5.012 X 10"'^ 


427.263 


Best fit point (approximate) 


Set 2 


0.031623 


1.995 X 10"^ 


428.806 


Large JIwdm good fit 


Set 3 


0.200000 


1.000 X 10-^ 


2425.883 


Sample point outside 2a contour 


Set 3* 


0.300000 


2.850 X 10"^ 


> 10000 


5^CDM = 0.004 



Table 1: [w, riwoM) points corresponding to figures (2), (3) and (4). 



power spectrum, both of which generate similar confidence regions. The supernova data places 
only rather weak constraints on our model, as it only probes low z where the model looks identical 
to ACDM. Our best fit model has a combined statistically equivalent to the x^ for ACDM. 
Figure (1) indicates the confidence regions in the w- ilwDM plane; figures (2), (3) and (4) depict, 
for the particular values of w and f^wDM listed in table (1), the redshift-distance modulus relation 
and the CMB power spectrum and matter power spectrum respectively. Unless otherwise noted the 
models assume standard values for cosmological paramters: f^wDE + ^wdm = 0.73, r^cDM = 
0.226 and il;, = 0.046. From these figures we see that there are parameter values for which smooth 
WDF provides as good a fit to all the available data as the concordance model based on ACDM. The 
trend from the x^ contours is to drive $7wdm to be relatively small, of order riwDM ^ 10^^ — 10" 
although larger values can be used by paying a relatively small price in the x? value. While the 
scaling behavior of WDM is such that it would dominate over radiation at some point, the smallness 
of r^wDM means that this event lies sometime in the future history of the Universe. 

What has become increasingly clear from various parmeterizations of the equation of state is 
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Figure 2: The distance modulus for a WDF dominated universe. 
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Figure 3: The CMB temperature anisotropy as computed for a WDF dominated universe. 



that it must be very close to ACDM at least up to z 1100. This follows primarily from the 
precision CMB and matter power spectrum measurements. It is important that the equation of state 
be physically reasonable and approximate ACDM in the appropriate regime. However, we need to 
be able to distinguish it from ACDM and hopefully this will be possible by considering the growth 
of perturbations. This assumes additional significance when we consider that the model differs 
significantly from ACDM only behind the surface of last scattering, a region quite inaccesible to 
current experiments. 



-4- 




Figure 4: The matter power spectrum in a WDF dominated universe. 



Let's now turn to the possibility that the WDF is also affected by the primordial density pertur- 
bations generated in the early universe, perhaps through an inflationary phase. If we consider the 
hydrodynamics of a fluid governed by some general equation of state characterized by p = w{p)p 
and v'^=dp/dp then linear perturbation theory tells us that the evolution of perturbations is governed 
by[29] 

5 + H5\2 - 3(2w - v"^)] - -H'^5[l - Qv'^ - Su-^ + 8it;] = - (—\ ^5. (10) 

2 \ a / 

where 6 = Sp/p and k is the wavenumber of the mode. In the ACDM case the cosmological 

constant is homogeneous so that we need only consider the perturbations in the CDM component 

for which eq.(lO) reduces to 

5 + 2H6 --H^6 = 0. 
2 

If we allow perturbations in the WDF, the scenario is quite different as the sound speed is a nonva- 
nishing constant and w{a) is no longer constant. At late times when w —1 we find 

6 + 8HS + 15H'^5 = -(^^y5. (11) 

However, the fact that w{a) does vary in time, will significantly modify the behavior of some of the 
modes. We recognize from the RHS of Eq. (1 1) a A; dependent suppression and that the coefficients 
of the homogeneous equations differ substantially from those of the ACDM case. We have used 
these equations in CMB-Fast to calculate the CMB angular correlations, and the matter power 
spectrum for various models. The results are quite interesting. If we normalize our calculated 
CMB power spectrum to the location and height of the first peak we find that even for ^^wdm ^ 
^^CDM, our model agrees well with the WMAP results as seen in figure (5). The change in \^ 
from the concordance model is just Ax^ = 5.3 for J^wdm = ^^cdm ■ Furthermore, with this 
normalization, our power spectrum also exhibits a low-/ suppression. 
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Figure 6: Matter power spectrum illustrating suppression. A tilted initial power spectrum can compensate 
for the suppression. 



As one might suspect, the matter power spectrum provides more restrictions on the model, 
as it probes higher values of k. As shown in Fig. (6) we observe an overall suppression int eh 
MPS, as well as an enhancement of the baryon oscillations leading to A^^ = 63.2. The matter 
power spectrum constrains ilcDM to be greater that r^wDM- However, this suppression can be 
compensated by increasing the tilt of the power spectrum without significantly affecting the fit to 
the WMAP data. This gives rise to a significant improvement of the fit, with Ax^ = 3.7. 
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We have considered both smooth and perturbed WDF fluids. However, it would be very useful 
to have a field theoretic, microphysical model that could give rise to the WDF equation of state. We 
could use this to study the evolution of perturbations in a more detailed manner and it's possible 
that entirely different properties may emerge when we consider candidate microphysical theories. 

Consider the following scalar field lagrangian [25]: 

C = C{X) where X = g^,u^^'(p^'' ^. (12) 

Such a Lagrangian can be made technically natural by the imposition of the shift symmetry (j) 
(j) + const. Computing the resulting stress-energy tensor we find 

T'^, = -C S^', + 2d''<l>d,cf>^, (13) 

from which we can obtain expressions for the energy density and pressure once we make assump- 
tions of spatial homogeneity and isotropy. These assumptions imply that di(f) = and X = <^^, so 
we can write 

p = C. (14) 

Now, assume that the "vacuum" state this theory finds itself in has a non-zero value for X[26]. 
Then we can arrive at a choice of £ that will give rise to the WDF equation of state by having C 
satisfy the following differential equation: 

2w ^^dC „ w 

-X— -C- = 0. (15) 



l + w dX l + w 

If we define 7 = and M an arbitrary constant with units of mass we can write the solution as: 

C{X) = {M^f-^X2 - ^. (16) 

7 

The Euler-Lagrange equation following from this Lagrangian is 

4> + 3w-^ = (17) 
a 

which impUes oc a~^'^. 

While this lagrangian does serve our purpose, further work is required to investigate the nature 
of perturbations about the homogenous solution. We begin this work by constructing the equation 
for field fluctuations around this background. This entails the replacing of (f) with (p + ^ and the 
evaluation of the contributions to the stress energy and equation of motion, keeping only the first 
order terms in ip. We find 

6T^ = -{j-l)% + p), (18) 

5Tf = -^{p + p) and (19) 
9 

ST^ = % + p)Si (20) 
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for the first order contributions to the stress energy, where the (p + p) are the background homo- 
geneous values. We need to choose a gauge to write the equation of motion and here we pick the 
conformal synchronous gauge[27] so that we find 

(7-l)v3+(4-7)-(^ + A;V + |# = 0, (21) 

where h represents the scalar metric perturbations. To go further we have to encode these contri- 
butions in CMBFast and integrate them. We will pursue this procedure in a future work. 

To conclude, we've started with an equation of state that is extremely simple and has a "mi- 
crophysical" parameter w that is positive and bounded so that the adiabatic sound speed in this 
fluid is less than 1. Furthermore, this equation of state describes the behavior of fluids as simple as 
water! From this we have extracted a cosmological model that agrees with all available data and 
is as statistically valid as the current favorite model, ACDM. We would argue that this ushers in 
a new paradigm in understanding dark energy in the sense that we can move away from demand- 
ing that the microphysical w be negative and we can certainly offer an alternative to considering 
w < —W Furthermore, we have come up with a proof of principle lagrangian that can describe this 
system.The next step is to see how the perturbations from this class of theories behaves, relative to 
ACDM, models, say. This work is in progress. 

Note: While we were writing this work up, [28] appeared which also deals with this equation 
of state. Their viewpoint is directed more towards the eventual fate of a Universe dominated by a 
WDF-like equation of state and does not have much overlap with the discussion in this paper. 
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